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Introduction 

Existen varios modelos que intentan definir la dinamica del desarrollo de tumores 
cancerigenos. Algunos involucran el crecimiento del tumor y la interaccion de celulas 
tumorales con virus oncoliticos (son virus que afectan y matan espedficamente celulas 
cancerigenas pero no normales). 

La viroterapia oncolitica puede facilitar la destruccion de las celulas tumorales cuando se 
administra quimioterapia (uso de drogas para tratar una enfermedad) y radioterapia. Es un 
tipo de terapia dirigida. Tambien se llama terapia virica, terapia virica oncolitica, y viroterapia. 
El efecto oncolitico resulta de al menos 3 tipos de interaccion huesped (ser atacado)-virus 
(agente causal) En cuyo caso los huespedes son las celulas cancerigenas. 

El primero involucra la replication ciclica viral en las celulas tumorales, que conlleva a la 
muerte de estas, consecuentemente a la reduccion del tumor y potencialmente a su 
eliminacion. El segundo involucra una reproduction de bajo nivel del virus que implica la 
produccion de la proteina citotoxina que causa dano a las celulas del tumor. 

El tercero involucra la infeccion de las celulas cancerigenas que causa inmunidad antitumoral. 
La infeccion viral causa inflamacion y penetracion de linfocitos al tumor y los antigenos 
virales (antigenos: cualquier sustancia que haga que el cuerpo parduzca anticuerpos) excitan la 
sensibilidad de las TNF para matar las celulas tumorales, abreviatura del ingles tumor necrosis 
factor, es una proteina del grupo de las citocinas secretadas por las celulas del sistema 
inmunologico y que tienen un efecto antitumoral 

A pesar que los ultimos 2 tipos de terapia son prometedoras, la primera es la estrategia mas 
efectiva actualmente, ya que no es indirecta, si no que se basa en lisis de las celulas cancerosas 
por un virus oncolitico. 

El factor de necrosis tumoral fue descubierto casualmente por William Coley (1893) que 
observo que en los pacientes con cancer que desarrollaban una infeccion bacteriana se inducia 
una necrosis tumoral. En la esperanza de curar el cancer, Coley trato a pacientes con cancer 
con el sobrenadante de algunos cultivos bacterianos. Las inyecciones de este sobrenadante, 
(llamadas toxinas de Coley) mostraron una cierta actividad en la necrosis tumoral, pero tenian 
efectos secundarios muy graves. Hubo que esperar hasta 1975 cuando se demostro que el 
sobrenadante de Coley tenia en realidad unas endotoxinas bacterianas que estimulaban a las 
celulas del sistema inmunologico a que secretaran una sustancia que si tenia un efecto 
antitumoral, siendo bautizada por este motivo como factor de necrosis tumoral. 

Wodarz describe la evolucion de tumores bajo inyeccion viral, presenta un modelo con 
interaccion entre celulas del sistema inmune y dos tipos de celulas tumorales (las que son 
infectadas por algun virus oncolitico y las que son susceptibles a la infeccion por ser 
cancerigenas). 
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El modelo siguiente esta basado en el modelo de Wodarz pero se considera solo la muerte 
directa de las celulas tumorales por efecto de un virus oncolitico sin tomar en cuenta el 
efecto celulas del sistema inmune. En contraste al modelo original de Wodarz, exhibe todas 
las posibles respuestas de la infeccion por virus oncoliticos. 


-^ = f 1(x ,y)X-8(X,Y)Y, 
clY 

-^ = f2 ( x,r)Y + g(X,Y)Y, 


f l{XJ) = r 1 (l-(X+Y)/K)-d, 
f 2 {XJ] =r 2 (l-(X+Y)/K)-a, 
g{X,Y) = bX 


Donde X(f) es la cantidad de celulas tumorales no infectadas y Y(f) la cantidad de celulas 
infectadas con el virus, r h r 2 , d, a, b, K son parametros no negativos. 

g(X,Y) es una respuesta funcional ,en este caso es el numero de presas consumidas por 
depredadores por unidad de tiempo para una determinada cantidad de presas X, 
depredadores Y , es decir, la cantidad de celulas no infectadas que son infectadas por unidad 
de tiempo por el virus que liberan las infectadas, similar a la ley de accion de las masas, se 
asume que la abundancia del virus es proporcional al numero de celulas infectadas y que el 
virus es heredado cuando una celula infectada se replica en 2,ya que los virus se integran al 
genoma de la celula, pero incluso si no, existe gran probabilidad de que se transmita a ambas. 

Se asume que fi(X, Y ) / = 1, es la tasa con la que crecen las celulas tumorales no infectadas por 
unidad de tiempo, la cual es una especie de funcion logistica donde d es la tasa de 
mortalidad natural de las celulas, r x es la tasa de crecimiento que se ve afectada por la 

capacidad de carga K del tumor, es decir, el numero maximo de celulas infectadas y no 
infectadas que pueden ocuparlo. 

Similarmente fi(X, Y) i =2 es la tasa con la que crecen las celulas tumorales infectadas por 
unidad de tiempo, la cual es tambien una especie de funcion logistica donde a es la tasa de 
mortalidad natural de las celulas, r 2 es la tasa de crecimiento que se ve afectada por la 
capacidad de carga K del tumor. 


Explicitamente el modelo es el siguiente: 


X- = {r\{\-{X + Y)/K)-cl)X -bXY, 

dY 

-X- = (r 2 (l-(X + Y)/K)-a)Y + bXY 



3 


Adimensionalizacion del modelo con d=0. 

La mortalidad natural explfcita d, puede ser eliminada sin perdida de generalidad (i.e., d = 0), 
ya que ambos sistemas tendran un retrato fase topologicamente equivalente, con la finalidad 
de tener el menor numero de parametros posible, con la misma finalidad tambien se puede 
escalar el tiempo y las variables que dependen de este. 

Dado d=0, sea: 1 dz = r x dt Entonces 

dX 1 

-= (1 -(X+Y)/K)X — bXY, 

dz r i 

— = -U(l-(X+Y)/K)-a)Y + - bXY 
dz r, n 

Sea: X — Kx dX — Kdx y Y = Ky dY = Kdy 


Kdx 

dz 

Kdy 

dz 


= (l-(x+ y))Ajr-— bK 2 xy, 
r i 

= -(r 1 (l-(X+Y))-a)Ky + - bK 2 xy 
r\ ' r i 


Y haciendo y = r-Jzi, 6 - bK/ri, 5 = a/ri. El sistema queda 


— = (1 ~(x+y))x~j3xy, 
dz 

X~ = (r( 1 -(x+y))-S)y + j3xy 


donde J3 > 0 


Para el cual se tienen los siguientes puntos fijos: 


7 ? =( 0 , 0 ) 

f 8\ 

P 2 = 0,1 — 

l YJ 

r, = ( i.o) 

p ( s+ps-Pr fl-s ' 
4 ypd+p-r)'pa+P-T)) 
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Los puntos de bifurcation que se obtienen son 
P 4 -P 3 = (1,0) cuando 5 = f3 


P, = P { = (0,0) cuando 5 = y 


P, 



1 ^ 
fi + ly 


cuando y = 8 


(yg + l) 

P 


s.q.^ * 13 


Retrato Fase- parametro del modelo adimensionalizado. 

El retrato fase- parametro se puede dividir en 6 dominios, cuyo analisis en principio 
sera sin contemplar sus fronteras, topologicamente diferentes del espacio parametrico 
(6, y, 6,) acotados por las fronteras ai, a 2 , a 2 , tal que: 


a 1 ={[6,y,6):6 = 6}, 


a 2 ={{6, y, 6): y = 6}, 


as={(6,y,6y.y = 6(6 + l)/6}. 
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Dado que para el dominio 1, 8 > P c\8 > y 


Dado que para el dominio II ,8 <y < — Or) > /? 


ns s $ $ i 

P <8 <y <8 -\ —n— >1 

/? P 


Dado que para el dominio III ,y >^^-^-r\8 > p 


e c e 

y >c> + —n£+ — >c>+l->y>c> + — >S+\ 

P P P 

~^y>S +1 
-^y-8 >l,s.q.S > P 


r> P 


Ademas 


Ya que y > 8+ \> 8 > P 


Dado que para el dominio IV, y > c\8< P 


8(J3 + 1 ) _ _ 


P 

y.Off+D 


— > 


/? 

X > SnS < p 


>1 nS < p 


Dado que para el dominio V, 8 <y r\y < ~ n $ < P 


s £ $ $ 1 

8 <y <8 + — n — < 1 
P p 



0 < y-8 < 1 ,s.q.8 < P 


Dado que para el dominio M\,y < 8 < P 
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Analisis de estabilidad del sistema 


Localmente a los puntos fijos el sistema se aproxima al siguiente 



= 3 

1 

X 

1 

X 

3 

_1 


2X2 

1 

e 

1 

_1 


, r dx c dy , \ 

Donde .£*= —, v*= —, ( x n , y n jes un punto fijo, 
dr dr 


3 = 


1 - 2-^ (l - v 0 - Py () 

-ry 0 +fiy 0 r~rx 0 -2yy 0 +^x 0 -s^ 


Para el punto fijo P l = (0,0) 



0 

y-8 


T rZs Pi =1 + y-S 

\z*\ = r-s 

Region I |3 P J<0 

Region II 13^, | > 0 

Region III 3,, | > 0 

Region IV |3 P |>0 

Region V |3 P J>0 

Region VI |3 P J<0 


Tr^ -4|3 / ,| = (l-r + ^) 2 >0 


TrS, 

Tr3, 

Tr3, 

Tr3, 

Tr3, 

Tr3, 


positivo o negativo pero no interesa. 

>0 

>0 

>0 

>0 

positivo o negativo pero no interesa. 
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Para el punto fijo P, = 0,1- 


i r-s i 
= o/ ,3 = 


£> + /?£> - Py 


ps 

P + r) -y8-y 




(S-r)(S+p(S- r )) 


T r3 p = 


(8-y)(p + y) + 8 


Tr %- 4 | 3 ,1= W±izll±r(-l±rl >„ 

2 I 2 


Como 


8>y ^ — > 1 -> y, < 0 


Por lo tanto este punto en la region I,VI es irrelevante. 


Region I 


3„ >0 Tr3 p >0 8 > y -> y 2 < 0 

r2 ' 2 ya que 


Region II 


3„ <0 Ya que (8 — y) <0,(8 + P(8 — y)) > 0 


Prueba: 


S(p + 1) i. 1 y 1 

y<- J - - >y<8 1 + — — <1 + — 

/? l * /? 


Tr3p Pero no interesa . 


Region III 3„ >0 Ya que (J -/) < 0 , (8 + P{8 - y)) < 0 


Prueba: 


8(P + 1) X 1) y , 1 

y>^~ — ~^y>8 1 + — — >1 + — 

/? l ^ /? 
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Tr3 p <0 T r3 p <-y D , 

^ ya que ' 2 ' Prueba: 


(S-y)(j3 + y) + S 

r 


< -y —> (S - y){j3 + y) + 8 < -y 2 —> (38 - py + yd + 8 < 0 


8{p + y + 1) - py < 0 —» S(J3 + y + 1) < J3y —> — (/? + ^ + 1) < y — > ^ > — (/? + 1) 


Region IV |3 P |>0 Ya que (<? - ;k) < 0 , (<5 + — y)) < 0 (la prueba se hizo con la 

region III). 


Tr3 D < 0 Tr3 D < — y (la prueba se hizo con la region III), 

ya que 


Region V 3^ | < 0 Ya que (8- y) < 0, (8 + [3(8- y)) > 0 
hizo con la region II). 

Tr3 

p * no interesa. 

Region VI |3 P |>0 Tr3 p >0 5 > y —> y 2 < 0 

I 2 1 12 ya que 


Para el punto P 3 = (l,0), 3 = 


-1 

0 


-1 -P 
P-5 


(la prueba se 


3 p =-j3 + S Tr3 p =-\ + j3-S 


Tr 2 3 p -4 


3 , 


= (-l-j3 + S) 2 >0 


En las regiones 1,11,111 donde S > J3, 3„ >0, Tr3 p < 0 


En las regiones IV,V,VI donde S < j3 ,3 P <0,Tr3„ puede ser positivo o negativo pero no 


interesa. 
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En la region I 
(J3-S)< 0 

135 + 5 - (3y > 0 Ya que 5 > y 
\ + (3 - y <0 ol + (3 -y >0 

Por tanto si 1 + (3- y < 0, x 4 < 0, y 4 > 0 en cuyo caso |3p 4 1 > 0, es decir P 4 e s un atractor. 
Similarmente si l + /3-y>0y 4 <0,x 4 >0 cuyo caso |3 P | < 0, es decir P 4 es un punto silla. 
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En la region II 

(P~8)< 0 

f38 + 5 - [3y >0 Prueba: 

y < 5 ^ 1} -> /3y<8(/3 +1) -> 0 < J35 + 5- J3y 
/?-X<0— >l + /?-x<0 ol + J3-y >0 

Por tanto si 1 + P~y < 0, x 4 < 0, y 4 > 0 en cuyo caso |3p 4 1 > 0, es decir P 4 e s un atractor. 
Similarmente si l + /3-y>0y 4 <0,x 4 >0 cuyo caso |3 P | < 0, es decir Z^es un punto silla. 

En la region III 

(J3~5)< 0 

J35 + 5 - (3y < 0 Prueba: 

r> (KP_p) -+p y> 8(p + \)^>Q > p8 + 8-py 

£ 

\ + P - y <Q Ya que \ + P~y<8 + — - f < 0 Prueba: 

a s 8 . 8/3 + 1 8(p + 1) 

— > 1 ^- — >—-> ——--> 6 +1 

P P P + l P 

c* c* 

-+8 + — >/? + l->£ + — -y > JS + l-y 

Por lo tanto x 4 > 0, y 4 > 0 |3p 4 1 < 0, es decir P 4 es un punto silla. 
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En la region IV 

iP~8)> 0 

P8 + 8 - Py < 0 La prueba se hizo en la region III. 

\ + P~y <0 o 14- P — y >0 

Por tanto si 1 + /?- y > 0,x 4 < 0, y 4 > 0 en cuyo caso |3p 4 1 < 0, es decir P 4 es un punto silla. 
Similarmente si 1 + P -y < 0 y 4 < 0,x 4 > Oen cuyo caso |3 P | > 0, es decir P 4 es un atractor. 

En la region V 

iP~8)> 0 

PS + 8 - Py > 0 La prueba se hizo en la region II. 

P(\ + P~ > 0 Ya que P(l + P — y) > P8 + 5 — Py Prueba: 

P> 8 P(P + 1) > 8(P + 1) p(p +1) > S(P +1) 

p 2 + P> PS + 5 —» p 2 + P~ Py > PS + 8—Py 

Por lo tanto x 4 > 0, y 4 > 0 |3p 4 1 > 0, es decir P 4 es un atractor. 

En la region VI 

iP~8)> 0 

PS + 8 - Py > 0 Ya que 8 > y 
P(l + P~y^>0 Ya que P>y 

Por lo tanto x 4 > 0, y 4 > 0 |3p 4 1 > 0, es decir P 4 es un atractor. 
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p ,=m F _ 

( 

0,1- — 

ii 

o 

ro 

ii 

(s+ps-fy P-5 ) 

1 2 

l rj 

4 

{Pd+P-rY PQ+P-r)} 


Region 1 

punto silla 

y<0,nodo 

fuente 

nodo 

estable 

x<0,y>0->nodo o espiral estables 
y<0,x>0->punto silla 

Region 

II 

nodo 

fuente 

punto silla 

nodo 

estable 

x<0,y>0->nodo o espiral estables 
y<0,x>0->punto silla 

Region 

III 

nodo 

fuente 

nodo estable 

nodo 

estable 

punto silla 

Region 

IV 

nodo 

fuente 

nodo estable 

punto silla 

x<0,y>0->punto silla 
y<0,x>0->nodo o espiral estables 

Region 

V 

nodo 

fuente 

punto silla 

punto silla 

nodo o espiral estables 

Region 

VI 

punto silla 

y<0,nodo 

fuente 

punto silla 

nodo o espiral estables 


La tabla resume el comportamiento del sistema en las diferentes regiones. 



Retrato Fase-parametro (tasa de mortalidad 0 para celulas no infectadas , 6 > 0, Y >0, 6>0). 

Las cotas entre dominios corresponden a cambios de comportamiento 
P 4 -P 3 - (l,0) cuando 5 = f3 Pr,-P x - (0,0) cuando 5 = y 

s.q. 5 * P 


P =P = 

r 4 1 2 


' 1 ^ 

0,- 

v P + b 


cuando y = 8 


(yg + D 

P 
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Analisis de estabilidad del sistema en los puntos de intersection de las cotas g-i, a? a 3 . 


Para el primer punto de intersection, i.e. ,5 = 15, 
resultante es: 


y = 5 


Qg + 1) 

p 


> y = (/3 + \) e\ sistema 


dx 2 .. .. 

- = x-x~ — (/? + l)xy, 
dr 

^ = - X y + y-(P + l)y 2 

dz 


Para el cual P, no existe. 


P 2 = 


y = 


f 1 > 

0,- 

v P + b 


yj p 3 =( 1 , 0 ) 


1 1 

-x + - 


/5 +1 J3 +1 


EnP 1= (0,0) se tiene una familia de rayos inestables, 
forman parte de un conjunto atractor dado por 


Suponiendo —— ,se sabe que el retrato fase sera una familia de rayos de la forma y = c l x, 

dy 

donde c x es constante, que pasan por el origen ya que el gradiente de 

dx x-x 2 -(B + l)xy 1 , 1 , 

— =--- es — ax - dy . 

dy -xy + y-(j3 + l)y x y 
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Para el segundo punto de intersection, i.e., S = J3 = / e I sistema resultante es: 


-y- = x-x 2 -(/3 + l)xy, 

Cl T 



Para el cual P 2 =/J=(0,0) se tiene una familia de Ifneas paralelas inestables, 
P 4 -P 3 = (l,0) se tiene una familia de Ifneas paralelas estables. 
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Conclusiones: 

Como se ha mencionado antes este modelo es un modelo simplificado del presentado por 
Wodarz , en tanto que ignora la dinamica de la poblacion del virus asf como la respuesta del 
sistema inmune, al tomarse en cuenta , se podrfa llegar a resultados mas confiables de la 
interaction del virus con el tumor. 

Si la tasa de mortalidad es nula para las celulas no infectadas... en la region I y II, 
eventualmente cualquier condicion inicial llegara al punto (K,0), lo cual quiere decir que no 
existiran celulas infectadas por el virus y se tendran K celulas no infectadas que es la capacidad 
maxima de carga del tumor, en otras palabras la terapia no tiene ningun efecto aunque se 
suministren altas dosis del virus. 

En la region III tambien se puede observar la misma situacion pero dependiendo de las 
condiciones infciales se podrfa llegar a un caso totalmente contrario, es decir donde no existen 
celulas no infectadas (todas estan infectadas) pero nunca se alcanzara la carga maxima del 
tumor, lo cual da la posibilidad de controlar los parametros para poder reducir el tamano del 
tumor. En la realidad no sucede que toda la poblacion de celulas se infecte por completo, esto 
puede ser debido mayormente a que no se consideraron algunos otros factores pero tambien 
se debe contemplar para que unidades de tiempo se tienen dichos parametros, pues es 
totalmente relativo el tiempo que debera pasar para que las poblaciones se vuelvan estables. 

Con el mismo criterio se puede analizar la region IV la cual se caracteriza por la 
predomination de celulas infectadas. 

En la region V y VI se puede observar que en general el sistema tiende a un punto de equilibrio 
donde tambien es posible reducir el tumor. Aunque en la realidad no se ha visto nunca el 
equilibrio entre ambas poblaciones, otra vez puede deberse a la falta de consideraciones (en 
los puntos de intersection entre las cotas tambien se puede observar un conjunto de puntos 
de equilibrio que dependeran de la condicion inicial). 

Los comportamientos en las Regiones I, II, III se deben entre otras cosas a que la tasa de 
infeccion es K veces menor a la tasa de muerte de las celulas infectadas por el virus, se puede 
asumir que las celulas infectadas mueren sin tener tiempo de infectar a otras, en las regiones 
IV, V, VI sucede que la tasa de infeccion es K veces mayor a la tasa de muerte de las celulas 
infectadas por el virus, por lo que se asume que las celulas infectadas infectan tantas antes de 
morir que eventualmente se llega a un equilibrio entre las dos poblaciones. 
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